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Vector dark energy and high-z massive clusters 
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ABSTRACT 

The detection of extremely massive clusters at z > 1 such as SPT-CL J0546-5345, SPT-CL 
J2106-5844, and XMMU J2235. 3-2557 has been considered by some authors as a challenge 
to the standard ACDM cosmology. In fact, assuming Gaussian initial conditions, the theoret- 
ical expectation of detecting such objects is as low as 1%. In this Letter we discuss the 
probability of the existence of such objects in the light of the Vector Dark Energy (VDE) 
paradigm, showing by means of a series of TV -body simulations that chances of detection are 
substantially enhanced in this non-standard framework. 
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1 INTRODUCTION 

Present day cosmology is still failing to explain satisfactorily the 
nature of dark energy, which is supposed to dominate the ener- 
getic content of the universe today and to be responsible for the 
current accelerated expansion. In the standard ACDM model, this 
cosmic acceleration is generated by the presence of a cosmological 
constant. However, the required value for that constant turns out 
to be tiny when compared to the natural scale of gravity, namely 
the Planck scale. Thus, the gravitational interaction would hence 
be described by two dimensional constants differing by many or- 
ders of magnitude, and this poses a problem of naturalness. This is 
the so-called "cosmological constant problem" and it motivated to 
consider alternative explanations for the current acceleration of the 
universe by either modifying the gravitational interaction at large 
distances or introducing a new dynamical field. 

Indeed, one of the main challenges of observational cosmol- 
ogy is exactly to devise new tests which could help discriminat- 
ing between the constant or dynamic nature of dark energy. In 
this regard, several authors have recently pointed out that the ob- 
servation of extremel y massive clusters at high redshift, such as 
SPT-CL 12106-5844 jFolev et alj bOllh . z ~ 1 18, M 20 o = 
(L27 ± 0.21) x 1O 15 M ), SPT-CL J0546-5346 terodwin et all 
J2010h . z ~ 1.07, M 2 po = (7 . 95 ± 0.92) x 1O 14 M ), and 
XMMU 12235.3-2557 jjee et alj | |2009|) . z ~ 1.4, M 200 = 
(7.3 ± 1.3) x 10 14 Mq) may represent a major shortcoming of the 
ACDM paradigm, where the presence of such objects should be in 
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princ i ple strongly disfavoure d (see, for example. lBaldi & Pettorinol 
1201 ltlMortonson et alj|201 ll) . 

While, on the one hand, this tension could be solved keep- 
ing the standard scenario and relaxing the assumption of Gaus sian- 
ity in the initial cond itions (as proposed in lHovle et alj fcOllh and 
lEnqvistetallfcOlllr ). it could be as well possible to use this ob- 
servations as a constraint for different cosmological models. In this 
work we look at the VDE model, wh ere the role of the dark en- 
ergy i s played by a cosmic vector field dBeltran limenez & Marotol 
2008). By means of a series of Af-body simulations, we study the 
large scale clustering properties of this cosmology, computing the 
cumulative halo mass functions at different redshifts and compar- 
ing them to the predictions of the standard model. In this way, we 
are able to show that the VDE cosmology does indeed predict a 
higher abundance of massive haloes at all redshifts, thus enhancing 
the probability of observing such objects with respect to ACDM. 



2 VECTOR DARK ENERGY 

The action of the vector dark energy model 
Beltran limene z~& Marotol J2008h ) can be written as: 



(see 
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and F^u = <9 M A„ 



g^ 1, Rfj,v the scalar curvature 
d v A^. This action can be interpreted as the 
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Maxwell term for a vector field supplemented with a gauge-fixing 
term and an effective mass provided by the Ricci tensor. It is in- 
teresting to note that the vector sector has no free parameters nor 
potential terms, being G the only dimensional constant of the the- 
ory. 

For a homogeneous and isotropic universe described by the 
flat Friedmann-Lemaitre-Robertson-Walker metric: 



ds 2 = dt 2 — a(t) 2 dx 2 



(2) 



we have — (Ao(t), 0, 0, 0) so that the corresponding equations 
read: 



A + 3HA -3 \2H Z +H\ A = 



2 SttG . , 

H = — — [p R + p M + Pa] 



(3) 



(4) 



(5) 



with H = a/a the Hubble parameter and 

p A = \a 2 Al + 3HA A - iig 

the energy density associated to the vector field, while pM and pn 
are the matter and radiation densities, respectively. During the ra- 
diation and matter eras in which the dark energy contribution was 
negligible, we can solve Eq. {3]l with H = p/t, where p = 1/2 for 
radiation and p = 2/3 for matter eras respectively, that is equiva- 
lent to assume that a oc t p . In that case, the general solution is: 



A (t) = A^t a + + Aot 



(6) 



with A constants of integration and a± = — (1 ± l)/4 in the 
radiation era, and a± = (—3 ± V33)/6 in the matter era. Af- 
ter dark energy starts dominating, the equation of state abruptly 
falls towards wde — > — oo as the Universe approaches t en d, 
and t he equation of state can cros s the s o-called phantom divide 
line jNesseris & Perivolaropoulosl l l2007t) ). so that we can have 
w D e(z = 0) < — 1. 

Using the growing mode solution in Eq. (JSJl we can obtain the 
evolution for the energy density as: 



PA = PAo(l + z) 



(7) 



with k = 4 in the radiation era and k — (9 — v33)/2 ~ —1.63 
in the matter era. Thus, the energy density of the vector field scales 
like radiation at early times so that the ratio pA / pn is constant dur- 
ing such a period. Moreover, the value of the vector field Ao during 
that era is also constant hence making the cosmological evolution 
insensitive to the time at which we impose the initial conditions 
(as long as they are set well inside the radiation dominated epoch). 
Also, the required constant values of such quantities in order to 
fit observations are pA/pi?| C ari y — 10 -6 and A c ^ rly ~ 10 _4 M P 
which can arise naturally during the early universe, for instance, as 
quantum fluctuations. Furthermore, they do not need the introduc- 
tion of any unnatural scale, thus, alleviating the naturalness or coin- 
cidence problem. On the other hand, when the Universe enters the 
era of matter domination, pA starts growing relative to pM eventu- 
ally overcoming it at some point so that the dark energy vector field 
becomes the dominant component. 

Once the present value of the Hubble parameter Ho and the 
constant A^" ly during radiation (which fixes the total amount of 
matter VLm) are specified, the model is completely determined. In 
other words, this model contains the same number of parameters as 
ACDM, i.e. the minimum number of parameters of any cosmologi- 
cal model with dark energy. Notice however, that in the VDE model 
the present value of the equation of state parameter too = —3.53 is 
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Figure 1. Equation of state of the vector dark energy model for the best fit 
to SNIa data, shown in the range of our simulations. 



radically different from that of a cosmological constant (cf. Fig.[T] 
where the redshift evolution of oj(z) is shown our the range of our 
simulations). Despite this fact, VDE is able to simultaneously fit 
supernovae and CMB data with comparable goodness to AC DM 
(Beltran Jimenez & Maroto 2008; Beltran Jimenez et al. 2009). 



3 THE DATA 
3.1 Simulations 

We wanted to estimate the probability of finding massive clusters at 
z > 1 in the VDE scenario compared to the ACDM one by means 
of CDM only A-body simulations. For this purpose, we chose to 
use a suitably m odified version of the publicly available GADGET- 
2 tree-PM code dSpringell2005) . which had to take into account the 
different expansion history that characterizes the two cosmologies. 
In Table[T]we show the cosmological parameters used in the differ- 
ent simulations. For the VDE model, we have used the value of Qm 
provided by the best fit to SNIa data; the remaining cosmological 
parameters have been obtained by a fit to the WMAP7 CMB data 
of the model, wo denotes the present value of the equation of state 
param eter of dark energy . For ACDM we used the Multidark Simu- 
lation dPrada et"ail 201 1 ) cosmologica l parameters with a WMAP7 
ag normalization i Larson et al .11201 lh . In addition, we also simu- 
lated a ACDM-vde model, which implements the VDE values for 
Qm and as in an otherwise standard ACDM picture. Although this 
model is certainly ruled out by current cosmological constraints, 
it provides nonetheless an interesting case study that allows us to 
disentangle the effects of these two cosmological parameters on 
structure formation in the VDE model. 

We chose to run a total of eight 512 3 particles simulations summa- 
rized in Tableland explained below: 



• a VDE (and a ACDM started with the same seed for the phases 
of the initial conditions) simulation in a 500 h' 1 Mpc box, 

• a second VDE (and again corresponding ACDM) simulation 
in a 1 ft _1 Gpc box, 

• two more VDE simulations with a different random seed, one 
in a 500 hT 1 Mpc and one in a 1 ft ~ 1 Gpc box, as a check for the 
influence of cosmic variance, 

• two ACDM-vde simulations in a 500 h" 1 Mpc and a 1000 
h' 1 Mpc box. 
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Table 1. Cosmological parameters for ACDM ACDM-vde and VDE. 



Model 






w 


CT8 


h 


ACDM 


0.27 


0.73 


-1 


0.8 


0.7 


ACDM-vde 


0.388 


0.612 


-1 


0.83 


0.7 


VDE 


0.388 


0.612 


-3.53 


0.83 


0.62 



Table 2. Af-body settings used for the GADGET-2 simulations, the two 
500/i — 1 Mpc and the two ih~ Gpc have the same initial random seed and 
starting redshift 2 s t a rt = 60 in order to allow for a direct comparison of 
the halo properties. The number of particles in each was fixed at 512 3 . The 
box size B is given in h~ 1 Mpc and the particle mass in /i _1 Mg. 



Simulation 


B 


nip 


2XVDE-0.5 


500 


1.00 x 10 11 


2 X VDE- 1 


1000 


8.02 x 10 11 


ACDM-0.5 


500 


6.95 x 10 1() 


ACDM-1 


1000 


5.55 x 10 11 


ACDM-0.5vde 


500 


1.00 x 10 11 


ACDM-lvde 


1000 


8.02 x 10 11 



The full set of simulations will be presented and analyzed in 
more detail in an upcoming companion paper; in this Letter, in- 
stead, we chose to focus on some of them only in order to gather 
information on large scale structures and massive cluster at high 
redshift, respectively, in the two cosmologies. To this extent, the 
use of the same initial seed for generating the initial conditions in 
the coupled ACDM vs. VDE simulations allows us to directly com- 
pare the structures identified by the halo finder. 

As a final remark, we underline here that the choice of the 
boxes was made in order to allow the study of clustering on larger 
scales, without particular emphasis on low mass objects, e.g. dark 
matter haloes with M < 10 14 /i _1 Mq. This means that even 
though our halo finder has been able to identify objects down to 
~ lO 12 /i _1 M in the 500/2T 1 Mpc box and ~ 10 13 /i _1 M in 
the 1 ft" 1 Gpc one (which corresponds to a lower limit of 20 parti- 
cles, see below), we are not comparing the mass spectrum at this 
far end. Therefore, since we are only interested in studying the be- 
haviour of the mass function of these models at the very high mass 
end, in the following sections we will mostly refer to the ACDM-1, 
ACDM-vde and VDE-1 simulations, where we have a larger statis- 
tics for the supercluster scales. 

3.2 Halo Finding 

In order to identify halos in our simulation we have run the 
MPI+OpenMP hybrid halo finder AHF described in detail in 
iKnollmann & Knebel J2009l> . AHF is an improvement of the MHF 
halo finder et alj |2004), which locates local overdensities in 
an adaptively smoothed density field as prospective halo centres. 
The local potential minima are computed for each of these density 
peaks and the gravitationally bound particles are determined. Only 
peaks with at least 20 bound particles are considered as haloes and 
retained for further analysis, even though here we focus on the most 
massive objects only. 

The mass of each halo is then computed via the equation 
M(r) = Ap c 47rr 3 /3 where we applied A = 200 as the over- 
density threshold. Using this relation, particular care has to be 



Table 3. The most massive halo found in the three \h 1 Gpc simulations 
(in units of 1O 14 /i -1 M0) as a function of redshift. 



z 


ACDM-1 


VDE-1 


ACDM-lvde 


1.4 


4.16 


5.63 


6.47 


1.2 


5.13 


6.51 


8.16 


1.1 


6.01 


7.63 


10.2 





18.1 


31.6 


35.1 



taken when considering the definition of the critical density p c — 
3H 2 /8tyG because it involves the Hubble parameter, that differs 
substantially at all redshifts in the two models. This means that, 
identifying the halo masses, we have to take into account the fact 
that the value of p c changes from ACDM to VDE. This has been 
incorporated into and taken care of in the latest version of AHF 
where Hvde(z) is being read in from a precomputed table. 

We would like to mention that we checked that the objects 
obtained by this (virial) definition are in fact in equilibrium. To this 
extent we studied the ratio between two times kinetic over potential 
energy 2T/\U\ confirming that at each redshift under investigation 
here this relation is equally well fulfilled for the ACDM and - more 
importantly - the VDE simulations (not presented here though). We 
therefore conclude that our adopted method to define halo mass 
in the VDE model leads to unbiased results and yields objects in 
equilibrium - as is the case for the ACDM haloes. 



4 THE RESULTS 
4.1 Mass Function 

With the halo catalogues at our disposal, we computed the cumu- 
lative mass functions n(> M) at various redshifts. We show in 
Fig.|2]the results for the lh, _1 Gpc simulations at redshifts z = 1.4, 
z = 1.2, z = 1.1 and z — 0. This plot is accompanied by Table[3] 
where we list the masses of the most massive haloes found in each 
model and the redshifts under consideration. 

We notice that the mass function for objecst with M > 
1O 14 /i _1 M is several times larger in VDE than in ACDM at all 
redshifts, i.e. the number of high mass haloes in this non-standard 
cosmological model is significantly increased. In particular, at this 
mass scale the VDE mass function is about three times larger at the 
relevant redshifts z = 1.4, 1.2, and 1.1 - and even larger at today's 
time. 

In order to verify that this feature of the VDE model is not a 
simple reflection of cosmic variance (which should affect in partic- 
ular the high mass end, where the statistics is small) we compared 
the results presented in Fig. [2] to the mass functions of the set of 
two additional simulations started from a different random seed for 
the initial conditions confirming aforementioned results. 

An interesting remark we would like to add here, is that the 
physical mass (obtained dividing by the corresponding h values the 
values quoted in /i _1 M units) of the largest haloes in the VDE-1 
simulation at z = 1.4, z = 1.2 and z = 1.1 are perfectly com- 
patible with the ones of the above clusters referred to in the Intro- 
duction, whereas the corresponding ACDM candidates are outside 
the 2a compatibility level. And again, similar massive clusters have 
also been found in the duplicate VDE-1 simulation with a different 
initial seed. 

Comparing the ACDM-lvde to the VDE-1 simulation at dif- 
ferent redshifts, we note that while the two mass functions are al- 
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most indistinguishable for M < 10 14 x /i _1 Mq, on the higher 
mass end the former even outnumbers the latter by approximately 
~ 3. In the hierarchical picture of structure formation, we can at- 
tribute this relative difference to dynamical effects caused by the 
different expansion histories (based upon different H[z)) at later 
times z fa 1, when the most massive structures actually start to 
form. In general, however, the ACDM-vde analysis shows that 
the enhancement seen in the VDE mass function with respect to 
ACDM is clearly driven by the higher values of Q.m and ag, a re- 
sult expected since the abundance of cluste rs sensitively depends o n 
the product of these two parameters (cf. iHuterer & Whitej[2002h . 
On the one hand, this complicates the issue of model selection, 
since (although disfavoured by the WMAP7 data) we could invoke 
a (slightly) larger Qm or a higher erg normalization at z = for 
ACDM trying to alleviate the current tension with the high-z mas- 
sive clusters observations. On the other hand, the distinct expan- 
sion history that characterizes and differentiates between the two 
ACDM and VDE models would still leave a clear imprint on struc- 
ture formation, which could be detected by, for instance, measur- 
ing as 's dependence on redshift. Such a test would indeed provide 
invaluable information for the study of ACDM and for any cosmo- 
logical model beyond it such as VDE. 

4.2 Probability 

In order to provide a more quantitative estimate of the the relative 
probability of observationally detecting such massive clusters at the 




0.5 1 1.5 2 

z 



Figure 3. Numerical cumulative number densities of objects with M > 
5 X 10 14 /l -1 M Q for VDE, ACDM and ACDM-vde. 



indicated redshifts we used n(> M, z) - the expected cumulative 
number density of objects above a threshold mass M as a function 
of redshift as given by our simulations - and integrated it over the 
comoving volume V c of the survey 

N(>M)= [ n(> M,z)dV c (z) (8) 

J Az,f2 survC y 
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Table 4. Expected number of objects JV(> M) in excess of mass M and 
inside a certain (comoving) volume in the ACDM and VDE for different 
mass thresholds and survey volumes. Solid angles Q are measured in deg 2 
and masses are measured in 10 14 /i — 1 Mq. 



M 


Az 


"^survey 


^Vacdm 


^VDE 




> 10 


> 1 


2500 


0.007 


0.02 


0.04 


> 7 


> 1 


2500 


0.03 


0.31 


0.56 


> 5 


1.38 - 2.2 


11 


0.005 


0.06 


0.07 



where Az and f2 sur vcy are the redshift interval and the fraction of 
the sky covered by the survey to which we want to compare our 
theoretical expectations. 

While n(> M, z) can be readily calcul ated in 

ACDM cosmologies (e.g. IPress & Schechted 1 1974 



ISheth & Tormenl 1 19991 : I Jenkins et al.l 1200 ll : iTinker et all [200 
in VDE we have to devise a strategy to compute it based upon 
our numerical results only. We chose to adjust the formula of 
ISheth & Tormenl i ll 9991) as follows: 

• we calculated the cumulative number densities in the desired 
redshift intervals Az based upon our simulation data, 

• we adjusted the parameters of the Sheth-Tormen mass func- 
tion fitting the numerical cumulative number densities derived from 
the VDE-1 and VDE-0.5 simulations, 

• we used these best fit estimates to analytically compute 
n(> M, z) now having access to masses even outside our numeri- 
cally limited range to be used with Eq. ^8). 

The results of the numerical integration over the comoving 
volumes ( obtained using t he limits quoted in t he o bservational 
papers bv ljee et al.l <2009b . iBrodwin et alj < f2oToh and iFolev et all 
( |20 111) ) are listed in Tablegjfor the VDE, ACDM-vde and ACDM 
model. We can clearly see that the chances are substantially larger 
in the VDE model to find such massive objects than in ACDM; the 
numbers for VDE are in fact comparable to our fiducial ACDM- 
vde model confirming their relation to the enhanced values of 
Qm and erg. However, note that while the VDE model is com- 
pliant with both SN-Ia and CMB data the ACDM-vde is obvi- 
ously ruled out observationally. We complement these results with 
Fig.[3]where we plot the abundance evolution of clusters with mass 
M > 5 x lO 14 /i _1 M0 computed utilizing above described pro- 
cedure again. This plot confirms our previous analysis of the mass 
functions and shows that the expectation of massive objects is am- 
plified in VDE by a factor ~ 3 to ~ 10 over the considered redshift 
range, a factor which is even higher for the discretionary ACDM- 
vde. We would like to remark here that while our ACDM estimate 
fo r XMMU J2235. 3-2557 is in agreement with the result quoted 
bv ljee et al. I d2009h (obtained using the same approach as here), the 
calculation done for SPT-CL J2 106-5 84 4 leads to an estim ate sub- 
stantially smaller than the one quoted bv lFolev et al.l d201 ll) . calcu- 
lated using a Monte Carlo technique. However, this does not affect 
our conclusions, which are based on the comparison of results ob- 
tained in a consistent manner for the two models. 



5 CONCLUSIONS 

The observation of massive clusters at z > 1 provides 
an additional, useful test for the ACDM and other cosmo- 
logical models beyond the standard paradigm. In this Letter 



we have shown that the Vector Dark Energy (VDE) scenario 
jBeltran Jimenez & Marotol 120081) might account for such obser- 
vations better than the ACDM concordance model, since the rel- 
ative abundance of extremely massive clusters with M > 5 x 
1O 14 /i -1 M0 is at all redshifts higher in this non-standard cos- 
mology: the expected number of massive clusters is enhanced in 
VDE by at least a factor of ~ 3 t o find an object s uch as SPT- 
CL J2106-5844 at redshift z « 1.2 jFolev et alj201lh and a factor 
of ~ 10 for the othe r two observed clusters SPT -CL J0546-534 6 
( IBrodwin et ai]|20ld) and XMMU J2235.3-2557 Jjee et al.ll2009b . 
Of course, these results might as well simply point in the direction 
of modifying the standard paradigm, for example including non- 
Gaussianities in the initial conditions or either using a higher erg 
or Qm value for the ACDM as the comparison to the ACDM-vde 
model seems to suggest. 

Nonetheless, this first results on the large scale clustering in 
the case of VDE cosmology point in the right direction, signifi- 
cantly enhancing the probability of producing extremely massive 
clusters at high redshift as recent observations seem to require. For 
a more elaborate discussion and comparison of the VDE to the 
ACDM model (not solely focusing on massive clusters) we though 
refer the reader to the companion paper (Carlesi et al., in prepara- 
tion). 
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